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Abstract 



We construct a new representation of composite fermion wave functions in 

the lowest Landau level which enables Monte Carlo computations at arbitrary 

filling factors for a fairly large number of composite fermions, thus clearing 

the way toward a more detailed quantitative investigation of the fractional 

quantum Hall effect. As an illustrative application, thermodynamic estimates 

for the transport gaps of several spin polarized incompressible states have 

been obtained. 
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Interacting electrons confined to two dimensions and subjected to a strong magnetic field 
exhibit spectacular phenomena, e.g., the fractional quantum Hall effect (FQHE) |I|]. A rather 
simple and succinct qualitative explanation of these phenomena is given in terms of objects 
called composite fermions 0, which are electrons bound to an even number of vortices of 
the many body wave function. The FQHE is a manifestation of the Landau level (LL) 
structure of composite fermions, and several recent experiments [H have produced striking 
additional evidence for composite fermions by detecting their semiclassical cyclotron orbits 
in the vicinity of the half filled LL [Q. 

A detailed quantitative description of the FQHE and related phenomena is less than 
satisfactory, however. It relies largely on exact diagonalization studies PJ^, which have 
played an extremely useful role in testing and confirming various theoretical postulates, 
but whose predictive power is rather limited as they typically deal with systems containing 
fewer than 10-12 electrons. Not only is this insufficient for making reliable thermodynamic 
estimates in most large region of filling factors is totally inaccessible in such small 

systems. To see this, take the conventional spherical geometry. Here, the n/{2n + 1) FQHE 
state requires at least electrons (which is the minimum number of composite fermions 
needed to fill n LLs). As a result, only two data points are available for 3/7 {N = 9 and 
12) in present exact diagonalization studies, and 4/9, 5/11, etc. do not show up at all. 
The situation is worse for the n/{4n + 1) sequence - even 3/13 is out of reach. Due to 
an exponential increase with the number of electrons in the computer time and memory 
requirements, it is unlikely that exact diagonalization studies will tell us much more in the 
future. 

A more promising approach toward the goal of a better quantitative description of the 
FQHE is to work with the composite fermion (CF) wave functions, which have been shown 
to be remarkably close to the exact solutions |^ and are expected to yield thermodynamic 
estimates for various quantities correct to within a few percent. Due to technical difficulties, 
however, it has not been possible in the past to work with these wave functions for more 
than ~ 10 composite fermions and, as a result, the CF theory has also in general failed 
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to provide quantitative information beyond what was known from exact diagonalization 
studies. In this Letter, we develop a new representation of the CF wave functions which 
enables computations for rather large CF systems. We report below results for as many as 
40 composite fermions; treatment of much bigger systems should be possible in the future. 
This constitutes a significant step in our ability to achieve detailed quantitative predictions 
for the FQHE. As a first application of this method, we have computed the transport gaps 
for several FQHE states of interest. 

The approach described below has been tested in both the disk (planar) and spherical 
geometries; here we will discuss only the latter due to space constraint. This geometry 



11,12] considers electrons on the surface of a sphere moving under the influence of a 



radial magnetic field. The single particle eigenstates are the monopole harmonics Yg^n 



TT| , which, after a bit of algebra, can be expressed as 



q,n,m 



(1) 



_,m / 2/ + l(/-m)!(/ + m)! \ 

where Qj represents the angular coordinates 6j and (pj of the jth electron, the total flux 
through the sphere is equal to 2g0o (2g is an integer and 0o = hc/e), n = 0, 1,2, ... is the 
LL index, m = — g — n, —q — n + 1, g + n labels the degenerate states in the nth LL, and 
I = q + n. It is understood here and below that the binomial coefficient (j^ vanishes if either 
/? > 7 or /5 < 0. The spinor coordinates are defined as [0 Uj = cos{9j /2) exp{—i(f)j /2) and 
Vj = sm{ej/2)exp{i(j)j/2) |13|. 

The CF theory postulates that the strongly correlated liquid of interacting electrons is 
equivalent to a weakly interacting gas of composite fermions. The {2p) vortices bound to 
electrons have the effect of partly cancelling the external field, and the composite fermions 
effectively experience a reduced magnetic field, given by -B* = -B — 2pp(j)Q, where B is 
the external field, and p is the electron density. Equivalently, the effective filling factor of 



composite fermions (z/*) is related to the electron filling factor u hj u = v* /{2pv* + 1) and, 
in particular, the FQHE of electrons at u = n/{2pn + 1) is a manifestation of the u* = n 
integer QHE (IQHE) of composite fermions. The (unprojected) wave function for the CF 
state at u*, ^cp, is related to that of the electron state at z/*, $, as 

= J'^, (3) 

where the Jastrow factor is given by 

J = YliujVk - VjUkf^exp[ip{(/)j + </)fc)] . (4) 
j<k 

One problem with is that it is, in general, not strictly in the lowest LL (LLL), and 
thus cannot be used as a good quantitative representation of the electron state at u in the 
high B limit. This is remedied by simply throwing away the part of the wave function that 
involves higher Landau levels and working with the remaining wave function V^Qp, where V 
is the LLL projection operator. It has been demonstrated to be very accurate by comparison 
with exact solutions available numerically for small systems 0. However, in spite of several 
attempts, it has not been possible to compute with V^'fjp for large CF systems, as would 
be necessary for exploiting the full quantitative potential of the CF theory. 

Let us outline the basic philosophy behind our approach. There are compelling theo- 
retical and experimental reasons to believe that ^cp contains the correct physics and is 
adiabatically connected to the true electron wave function in the lowest LL. Moreover, it 
has only a small fraction of electrons outside of the lowest LL Therefore, the objective 



is to obtain a LLL wave function from ^^p without disturbing it violently. There is noth- 
ing a priori to choose from between various methods that accomplish this goal; of course 
the accuracy of the resulting LLL wave function must be established by comparison with 
exact solutions. The principal result of this work is that we have constructed a LLL wave 
function for composite fermions different from V^Qp but comparably accurate and, most 
importantly, much easier to compute with. 

For simplicity, we will consider below only states for which $ is a single Slater determi- 
nant, $ = Det[Yi{Qj)]; generalization to the case where $ is a linear superposition of Slater 



determinants is straightforward. Then, 

^'cF = J Det%{n,)] . (5) 

We write the Jastrow factor as 

^ = n ("i^fc ~ VjUkY exp[i-(0j + 0fc)] = n ' (6) 

with 

■^j = n("i^fc ~" ^i^fc) exp[-((/'j + (t)k)] , (7) 
fe ^ 

where the prime denotes the condition j ^ k. Jj has the property that when expanded in 
terms of a hnear superposition of single particle eigenstates of the jth electron, all eigenstates 
correspond to the same monopole strength q' = {N — l)/2 (although with different m'). The 
Jastrow factor can now be incorporated into the Slater determinant to give 

%F = Det[Y,{Q,)Jj;] . (8) 

Now, instead of projecting the determinant on to the lowest LL, as had been done previously, 
we project each matrix element individually to write 

= Det[PY,{n,)J^] . (9) 

Since the product of two LLL wave functions is also in the lowest LL, ^cf is guaranteed to 
be in the lowest LL. In order to evaluate the projection VYq^n^miS^j) Jj, we first show that 
there exists an operator Y^'„ ,^ satisfying the property that 

Vyq,n,myq',0,m' = Y^,n,m^g',0,m,' , (10) 

where Vg'^cm' ~ e^'^''^^Uj '^'^ vj is a LLL wave function at monopole strength q'. For the 
present purposes, it is important that Y^'„ be independent of m'. To this end, we multiply 
one of the terms on the right hand side of Eq. ([I|) by the LLL wave function Fg'.o,m' and 
write (with Q = q + q',M = m + m', and the subscript j suppressed): 



e'^'^{v*v)''-'{u*uyu^+^v^-^ = aoe'^^'^u^+^v^-^ + higher LL states. (11) 

For \M\ > Q, ao must vanish, since \M\ < Q in the lowest LL. Let us first consider the 
case \M\ < Q. Multiphng both sides by Q-'^Q<t'ii*Q+^^y*Q~'^^ and integrating over the angular 
coordinates gives 

^ {Q-M + n- s)\{Q + M + s)\{2Q + 1)! 
"° {Q + M)\{Q-M)\{2Q + n + l)\ ' ^ ^ 

This shows that, apart from an m'-independent multiplicative constant {2Q + 1)1/ {2Q + 
n + 1)!, the LLL projection of the left hand side of Eq. (|Tl|) can be accomplished by first 
bringing all u* and v* to the left and then making the replacement 

^ du ' ^ dv ' ^^'^^ 
While this prescription is not valid in general for \M\ > Q, it can be shown to produce 
the correct result (i.e., zero) even for \M\ > Q for the LLL projection of states of the form 
^g,n,m^q',o,m': tcrms with uou-zcro binomial coefficients have the form (^^^ (^-^^ 
with either a > (3 (for M < Q) or 7 > 5 (for M > Q), and consequently vanish. Thus, 

V?' — AT (_T\q+n~mr)rn iq(j)j 



V 



n—s 



A delightful simplification occurs when one brings all the derivatives to the right in Eq. (|T4D 
using 

and a similar equation for the derivative with respect to u (with the summation index a'). 
The sum over s in Eq. (^4|) then takes the form 

i:Vi)-(":!7%"r'(-ir-(""r1- ^^^^ 

s=a' \ / s'=0 \ / 

which is equal to (— 1)"'(1 — i)"~°-°' and vanishes unless n = a + a'. The only term 
satisfying this condition is one with a = n — s and a' = s. Consequently, the derivatives 
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in Eq. ([M]) can be moved to the extreme right and act only on the following LLL wave 
function. Specializing to the case where Y^'„ ,„ acts on Jj, we write 

' j^=j^Rr-^ , (17) 



^ duj J \ dvj 



where 



= Jr^Jj = + 7^ ' (19) 

Y . = J7P—JP = pY h— . (20) 

Substituting in Eq. (P) and factoring out the mini-Jastrow factors Jj to produce back the 
full Jastrow factor J' finally gives 

<!>CF = J Det[Y,{Q,)] , (21) 

which has the same form as the unprojected wave function ^'cp in Eq. (^) except that each 
Y has been replaced by Y, given by 

Y (OA — N ( T\g+n-mom (2Q + 1)! q+m q-m 

|:(-i)'(:)(,,':::_>rs-«r'. 

A consequence of the LLL projection is that changing the coordinates of one particle alters 
all elements of the Slater determinant in $cf, as Rj depends on all particle coordinates. As 
a result, certain time saving tricks for updating the Slater determinant |]15| at each step of 



the Monte Carlo cannot be used here and the full determinant must be evaluated at each 
step. This increases the computation time enormously (as compared to Monte Carlo on 
^'cf)j but it is still possible to deal with much bigger systems than before. 

Of course, the usefulness of ^cf hinges critically on how precisely it approximates the 
true electron state. We resort to small systems to answer this question. Here and below, we 



will assume fully polarized electrons confined to the lowest LL. The Coulomb energy of $cf, 
obtained by Monte Carlo, along with the corresponding exact ground state energy is given 
in Table I for several incompressible states. A comparison establishes the extreme accuracy 

of $CF- 

We are now in a position to make detailed quantitative predictions for a number of 
experimentally measurable quantities in the FQHE. We consider in this paper the (transport) 
gaps of FQHE states belonging to the principal sequence v = n/{2n + 1). The gaps of the 
1/3 and 2/5 states have been known quite well from exact diagonalization studies [^^; the 
gap of 3/7 is known with a large (40%) uncertainty [§; and no estimates exist for 4/9, 5/11, 
etc. From a different perspective, Halperin et al. 0] have suggested, based on particle-hole 
symmetry in the lowest LL and the earlier known numerical results, that the gaps of the 
n/ {2n + 1) states are given by the equation 

C 

E„= , , ^ , (23) 

^ 2n + l ei ^ ^ 



with C ~ 0.31. There is also experimental support for such behavior ||T6[|. Equating Eg to 
the effective cyclotron energy of composite fermions, HeB* /■m*c, motivated by a mean-field 
description of composite fermions, produces an effective mass of composite fermions that 
scales as m* ~ ^/B [|]. 

The gap of the FQHE state at z/ = n/{2n + 1) is equal to the energy required to create a 
far separated CF particle-hole pair. The ground state has n filled LL's of composite fermions 
and the excited state is obtained by taking one composite fermion out from the south pole 
of the nth CF-LL and placing it on the north pole of the {n + l)th CF-LL (to maximize 
the distance between the CF-particle-hole pair). Since the gap, an 0(1) energy, is obtained 
as the difference between two large 0{N) energies, its accuracy is expected to be much less 
than that of either the ground state or the excited state energy. To get a feel for how well the 
CF theory does here. Table 11 quotes the gaps predicted by the CF theory for some of the 
largest systems for which exact gaps are also known; a comparison indicates that the error 
in the thermodynamic limit is expected to be within a few percent, typically much smaller 
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than the statistical uncertainty of our Monte Carlo. We proceed to compute the gaps for 



2/5, 3/7, 4/9 and 5/11 as a function of |T^. Fig. 1 shows the results for 3/7 and 4/9. 
Each error bar signifies one standard deviation, as determined from the average gaps from 
ten Monte Carlo runs. The thermodynamic estimates, shown in Fig. 2, have been obtained 
using a chi-square fitting that biases the points by their error bars. For z/ = 3/7 the gap is 
consistent with the earlier estimate, although with much smaller uncertainty. The overall 
behavior of the gaps is in good agreement with Eq. (p3D , as seen in Fig. 2. It should be noted 
that the experimental values of the gaps are reduced by various unavoidable effects, e.g.. 



finite width of the quantum well [|l^], LL mixing ||T^ and disorder, and the present estimates 
are actually only the upper bounds. The finite well width softens the Coulomb interaction 
at short distances, which is straightforward to incorporate into our Monte Carlo and will be 
the subject of a future work. LL mixing is harder to deal with. Fixed phase Monte Carlo 
pO| or a variational approach considering a linear combination of ^'cf ^cf 0| may be 
useful in this context. 

In summary, we have developed an approach that permits an investigation of large CF 
systems and thereby makes accessible previously unexplored regions of the FQHE. We believe 
that it will prove extremely useful toward a better quantitative understanding of the FQHE. 

This work was supported in part by the National Science Foundation under Grant no. 
DMR93-18739 and by a fellowship from the John Simon Guggenheim Foundation. 
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Table Caption 

Table 1. Comparison between the energy per particle of the CF wave function ^cf and 
the exact Coulomb energy for the ground states at z/ = 2/5 and 3/7. The energies are in 
units of e^/ei, where i is the magnetic length and e is the background dielectric constant, 
and include interaction with the uniformly charged positive background. The energy of ^cf 
has been evaluated by Monte Carlo, with the statistical uncertainty in the last two digits 
shown in brackets. The exact results for = 12 and u = 3/7 are taken here and in Table 
11 from S. He, S.H. Simon and B.l. Halperin, Phys. Rev. B 50, 1823 (1994). 

Table 11. The gap to create a CF-particle-hole excitation in which a composite fermion 
is removed from the south pole of the topmost filled CF-LL and placed on the north pole of 
the lowest empty CF-LL. The exact gaps are also given; the gap at u = 2/5 for = 10 has 
been read off from the spectrum in Ref. [^]. 

Figure Captions 

Fig. 1. The energy gaps (in units of e^/ei) at 3/7 and 4/9 for several A^. In order to 
minimize finite size effects, the interaction energy of two appropriately charged pointlike 
particles situated at the two poles has been subtracted from the gaps here (compared to 
those quoted in Tables 1 and 11), following 

Fig. 2. Thermodynamic values of the gaps plotted as a function of l/{2n + 1). The 
gap for u = 1/3 is taken from Ref. 0]. The straight line is a chi-square fit, given by 
Eg = -0.001(5) + 0.320(2)(2n + 1)-^. 
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V 


N 


CF energy 


exact energy 


2 

5 


6 


-0.500339(42) 


-0.5004002 




8 


-0.480216(33) 


-0.4802436 


3 
7 


9 


-0.499138(71) 


-0.4991843 




12 


-0.482507(49) 


-0.4826388 



TABLE I 



V 


N 


CF gap 


exact gap 


2 

5 


6 


0.07615(61) 


0.07505 




8 


0.07021(87) 


0.06809 




10 


0.0681(12) 


0.0673(7) 


3 
7 


9 


0.0691(14) 


0.0681 




12 


0.0518(12) 


0.0525 



TABLE II 
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